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Abs tract. We present an improved version of the method of photometric mode identification of Heynderickx et 
al. (1994). Our new version is based on the inclusion of precise non-adiabatic eigenfuncti ons de termined in the 
outer stellar atmosphere according to the formalism recently proposed by Dupret et al. ( [2002a ) . Our improved 
photometric mode identification technique is therefore no longer dependent on ad hoc parameters for the non- 
adiabatic effects. It contains the complete physical conditions of the outer atmosphere of the star, provided that 
rotation does not play a key role. We apply our improved method to the two slowly pulsating B stars HD 74560 
and HD 138764 and to the /3 Cephei star EN (16) Lac. Besides identifying the degree £ of the pulsating stars, our 
method is also a tool for improving the knowledge of stellar interiors and atmospheres, by imposing constraints 
on parameters such as the metallicity and the mixing-length parameter a (a procedure we label non-adiabatic 
asteroseismology) . 
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1. Introduction 



Stars: atmospheres - Stars: variables: /3 Cep - Stars: variables: Slowly Pulsating 



A crucial problem in asteroseismology is mode identifica- 
tion. Firstly because, from a theoretical point of view, de- 
spite the linear non-adiabatic predictions, the mode selec- 
tion mechanisms are not well understood for many kinds 
of pulsating stars (S Scuti, (5 Cephei, slowly pulsating 
B stars, 7 Doradus, roAp stars, ...). Secondly because, 
from an observational point of view, we do not resolve 
the disks of stars other than the Sun so that we can only 
observe disk-integrated quantities. Thirdly, the rotational 
splittings and the "avoided crossing" effect produce such 
a complicated power spectrum that a mode identification 
based on the frequencies alone is generally impossible. 

Currently, both spectroscopic and photometric mode 
identification techniques are being used. The latter meth- 
ods are based on multi-colour photometry, and are the 
subject of this paper. The principle of these methods is 
to observe the photometric variations due to stellar os- 
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cillations in different colours and compare them to the 
theoretical predictions at the appropriate wavelengths 
(Dziembowski 1977b), St amford & Watson 
1988 L Garrido et al. |1990| , Heynderickx et al. 
2000). However, all these methods have 
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Watson 
1994 Garrido 
an important 

drawback: their theoretical predictions are very sensitive 
to the non-adiabatic temperature variations at the photo- 



sphere (Cugier et al. 1994, Balona & Evers 1999) 



Dupret et al. (2002a) developed a non-adiabatic code 
including a detailed treatment of the pulsation in the outer 
atmosphere. In this paper, we show how this non-adiabatic 
treatment opens the way to a significant improvement 
of the discriminant power of mode identification meth- 
ods based on multi-colour photometry. Indeed, by using 
Dupret et al.'s calculations, we are able to eliminate the 
weakest point of the mode identification method, namely 
the ad hoc parameters to express the non-adiabaticity. We 
present the results obtained for two slowl y pul sating B 
stars (SPBs) observed by De Cat & Aerts ( |2002|) and for 
the P Cephei star EN (16) Lac. 



2. Monochromatic magnitude variations of a 
non-radially pulsating star 

Theoretical expressions for the monochromatic magnitude 
variations of a non-radially pulsating star have been de- 



2 



M.-A Dupret et al.: An improved method of photometric mode identification 



rived by different authors. Dziembowski (1977b) was the 
first to derive an expression for the bolometric magni- 
tude variation of a non-radially pulsating star. He sug- 
gested also that a Wesselink technique could be formu- 
lated from these expressions. Balona & Stobie (1979) re- 
cast the suggestion of Dziembowski in an observationally 
feasible way. Stamford & Watson ( 198l| ) derived an ex- 
pression for the monochromatic magnitude variations of 
a non-radially pulsating star. They proposed to compute 
the local emergent monochromatic flux variation on the 
base of equilibrium atmosphere models (see Eq. (|J) be- 
low) and they simplified the way to compute the influence 
of the stellar surface distortion. Watson (1988) improved 
the method by taking the variation of the limb darken- 
ing into account (see Eq. (||) below), and discussed the 
importance of the different terms in the equation giving 
the monochromatic magnitude variat ions of a non-ra dially 
pulsating star. Garrido et al. ( 199C ) and Garrido ( 2000| ) 
derived a method of mode identification using Stromgren 
photometry, based on the formalism of Watson ( 1988 ) , and 
applied it to 6 Scuti and 7 Doradus stars. Heynderickx et 
al. ( |1994[ ) derived an expression for the surface distortion 
of a non-radially pulsating star in a Lagrangian formalism. 
He developed a method of mode identification based on 
photometric amplitude ratios and applied it to (5 Cephei 
stars. 

In all the previously cited papers, the non-adiabatic 
character of the pulsation was neglected or treated with 
an ad hoc parameter. Cugier et al. ( |1994[ ) were the first 
to use non-adiabatic computations for photometric mode 
identification, using the non-adiabatic pulsation code of 
Dziembowski ( 1977a ). The same code was also used by 
Balona & Evers ( 1999 ) for mode identifications of 5 Scuti 
stars. Townsend ( 2002 ) used a non-adiabatic code for the 
photometric modelling of SPBs. 

It is useful to detail the assumptions made, sometimes 
implicitly, by the previous authors and in our method, 
in order to derive an expression for the monochromatic 
magnitude variations of a non-radially pulsating star. In 
sections [| and ^, we will put forward the improvement 
and specificities of our method, compared to the one of 
the previous authors. 

al) We work in the linear approximation. 

a2) We neglect the coupling of modes due to the inter- 
action between rotation and pulsation. The angular 
dependence of a non-radial mode is thus described by 
a single spherical harmonic. 

a3) We assume that the gas column of the atmosphere 
at a given angular position (0, <f) is well described by 
a plane parallel atmosphere, which we call the local 
atmosphere. 

a4) For the geometrical distortion of the stellar surface, 
we work in the one-layer approximation. It is assumed 
that the visible part of the star, i.e. the photosphere, 
can be described by a single surface which is spherical 
at equilibrium. The radius Rq of this sphere is the ra- 
dius of the star, and in our method, we assume that it 



corresponds to the layer where the local temperature is 
equal to the effective temperature of the star. During 
the pulsation, it is assumed that this surface follows 
the movement of the matter. The surface distortion 
can therefore be deduced from the displacement field: 

a5) We assume that, during the pulsation cycle, the 
monochromatic outwards flux F\ + of the local at- 
mosphere is, for each given time, the same as the 
monochromatic outwards flux of an equilibrium plane 
parallel atmosphere model. 

a6) We assume that, during the pulsation cycle, F\ + re- 
mains perpendicular to the photosphere. 

a7) We assume that F\ + does not depend on the optical 
depth in the local atmosphere. 

a8) We assume that the local atmosphere depends only 
on two varying parameters: the local effective temper- 
ature and the local effective gravity. The chemical com- 
position of the local atmosphere is assumed to remain 
constant. 

a9) We assume that during the pulsation cycle, the limb 
darkening law h\ of the local atmosphere is, for each 
given time, the same as the limb darkening law of an 
equilibrium plane parallel atmosphere model with the 
orientation given by assumption (a4) and (a6). 

Assumption (a4) needs some precisions. We have 
shown in Dupret et al. ( 2002a ) that the relative 



Lagrangian variation of the optical depth (5t/t) is not 
negligible in stellar atmospheres. However, on the basis of 
non-adiabatic computations, we have checked that the rel- 
ative difference between the displacement of constant opti- 
cal depth layers and the "real" displacement of the matter 
(St/(kpRq)) is very small. For example, \6t\/(k p Rq) ~ 
0.005 at the photosphere of a typical [3 Cephei model and 
for the fundamental radial mode, with a relative radial dis- 
placement normalized to 1 at the photosphere. Therefore 
assumption (a4) is appropriate for the determination of 
the geometrical distorsion of the stellar surface (at least 
for g-modes and moderate order p-modes). However, we 
showed in Dupret et al. (2002a) that it is not appropriate 



to assume that the Lagrangian variation of the tempera- 
ture is equal to the variation of the temperature at con- 
stant optical depth, because dlnT/dr \St\ is not negligible 
compared to \ST/T\ at the photoshere. 

We note that assumption (a7) concerns only the flux. 
We do not make this hypothesis for the temperature which 
depends strongly on the optical depth in stellar atmo- 
spheres (see Dupret et al. [2002a , Fig. 1). 

Assumptions (a5) and (a9) have the same physical jus- 
tification as explained in Dupret et al. (2002a). Because of 
the very small thermal relaxation time of the atmosphere, 
we assume that, at each time of the pulsation cycle, the 
local atmosphere remains in radiative equilibrium. Using 



this assumption, Dupret et al. (2002a) assumed that, dur- 
ing the pulsation cycle, the temperature distribution (T(r) 
law) in the local atmosphere was, for each given time, the 
same as the temperature distribution of an equilibrium at- 
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mosphere model. For the same physical reasons, we make 
now the same assumption for the local monochromatic 
outwards flux F\ + and for the limb darkening law h\. 

From assumptions (a5), (a7) and (a8), the monochro- 
matic outward flux in the local atmosphere is given by: 



F 



-) + 6F+(9,<f>,t) 



F+ [ (T cff ) + ST cff (9, </>, t), 5o + Sg e (9, 0, t) ] ,(1) 



where = 



\F\ + \. In the linear approximation, we have 
thus by perturbing Eq. (|l]): 



dlnF{ 
d\nT eS 
ST eS 



T, 



off 



ST cS 
T fr 



d\nF+ 
ding 



$9c 
9c 



(2) 



Eq. (|I|) was first proposed by Stamford & Watson ( 1981j ). 
We proceed similarly for the variation of the monochro- 
matic limb darkening. From assumptions (a4), (a6), (a8) 
and (a9), we obtain in the linear approximation: 



hx 



d\n h\ 
ainr eff 
d\nh\ 



cir 



T, 



cir 



d In h\ 
d\ng 



Sg c 

9c 



(3) 
(4) 



where h\ is the normalized limb darkening law 

(Jo ^(A 4 ) A* d/i = 1), /x = ro • e z /, n is the normal to 
the photosphere and e z / is the unit vector pointing to- 
wards the observer. A similar equation was first proposed 
by Watson p88| ). 

We denote by e the amplitude of relative radial dis- 
placement in the photosphere: 



£ r (6,<j>,t) = R ePP (cos 9) cos(at + m4>) , 



(5) 



where 9 and <fi are the usual spherical coordinates with re- 
spect to the polar (rotation) axis of the star, P™(x) is the 
associated Legendre function of degree £ and azimuthal 
number m and er is the angular oscillation frequency. 
The quantities ST e s/T e s and 5g e /g c can be com- 



puted by the non-adiabatic code of Dupret et al. ( |2002a| ). 
Because the Eulcrian variation of the gravitational poten- 
tial at the surface is always very small ( <C \g£r\ ), 
it appears that in very good approximation, Sg c /g c is in 
opposite phase with the radial displacement. These two 
quantities can be expressed in term of the associated 
Legendre functions: 



cir 



T e ff 

6gc 
9c 



(0,<j),t) = f T ePp (cos 9) cos(at + m<f> + ip T ) 
(0,<p,t) = - f g ePp (cos 9) cos(at + m<p) , 



(G) 

(7) 
and 



where fx and f g are the amplitudes of ST^/T^ 
Sg c /g c corresponding to a normalized radial displacement 
at the photosphere. On the basis of the previous assump- 
tions and equations, an expression can be derived for the 
monochromatic magnitude variation of a non-radially pul- 
sating star: 



2.5 

5m\ = - — — e PP(cosi) bn 
In 10 

- (£-l)(£ + 2) cos(cri) 
+ fr cos(at + ip T ) (a T \ + Pt\ ) 

~ fg COs(tJi) (a g \ + /? S A ) ) , 



(8) 



where 5m\ is the variation of the monochromatic mag- 
nitude as seen by the observer, i is the inclination angle 
between the stellar axis and the direction towards the ob- 
server and 



(3tx 



hxfiPedfj,, 
o 

d In bex 
d\nT cS ' 
d In bjx 
ding 



(9) 
(10) 
(11) 



For the derivati ons le ading to Eq. (g), we refe r 
to Dziembowski ( |l977bj ), Stamford fc Wa tson ( |l98l[) , 



Watson ( |1988|) and Heynderickx et al. (|1994|) . 

In Eq. (^), the term proportional to (£ — 1)(£ + 2) cor- 
responds to the influence of the stellar surface distortion, 
the term proportional to fx corresponds to the influence of 
the local effective temperature variation and the term pro- 
portional to f g corresponds to the influence of the effective 
gravity variation. In our applications, we computed the co- 
efficients arx and a g x (derivativ es of t he monochromatic 
flux) from the models of Kurucz Jl993 ). An analytical law 
for the limb darkening is needed for the computation of 
bix and its derivatives. For the pr esent paper, we used 
a quadratic law (Wade & Rucinski 1985 ). We note that 
an improved non-linear limb darkening law has been pro- 
posed by Claret ( |2000| ) , but his computations were only 
done for Stromgren filters, while our applications concern 
Geneva and Johnson filters. 

In multi-colour photometry, one observes the integral 
of the monochromatic magnitude variation over the re- 
sponse of the filter: 



<5m, = 



XCT W *( A ) dA 



dA 



(12) 



where Wi(X) is the response curve of the filter i. Therefore, 
the different terms of Eq. (|8) depending on A have to be 
integrated, following Eq. (|12). 

3. Our version of the mode identification method 

The linear theory does not permit to predict the ampli- 
tudes of the eigenfunctions. Therefore, it is appropriate to 
use amplitude ratios and phase differences between differ- 
ent filters when comparing the theoretical predictions to 
the observations. On one hand, for S Scuti stars, the ob- 
servations and the theoretical predictions give significant 
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phase-lags between the different filters. Mode identifica- 
tion methods usin g the se phase-lags have bee n prop osed 
by Garrido et al. ( |1990D and Balona & Evers fll999| ). On 
the other hand, for (3 Cephei stars, SPBs and 7 Doradus 
stars, no phase-lags are observed (in agreement with the 
small phase-lags predicted by the theory). For the latter 
stars, mode identification methods based on amplitude ra- 



tios are thus appropriate (Heynderickx et al. 1994). We 
adopt here such a method. The theoretical procedure of 
our mode identification method is the following: 

1. We compute a stellar model with the appropriate ef- 
fective temperature, luminosity and mass. In our ap- 
plications, we used the new Code Liegeois d'Evolution 
Stellaire written by one of us (R. Scuflaire). 

2. We perform non-adiabatic computations for different 
degrees £ and for pulsation frequencies close to the 
observed ones. In our applicati ons, w e used the non- 
adiabatic code by Dupret et al. ( 2002a ). These compu- 
tations give the coefficients Jt, V't and f g for different 
degrees £. 

3. For each filter j and for each £, we compute 



.4 



j,th 



(1 -£)(£ + 2) 
-f T e^ T (or j + 0Tj) 
' fg + Pgj) 



(13) 



using the atmosphere models of Kurucz (1993) to com- 
pute o>t and a g . 

We choose a reference filter (indicated with subindex 
1). For B stars, this reference filter is the U filter 
giving the highest amplitudes and thus the highest 
S/N ratio. We compare the theoretical amplitude ra- 
tios (A,,thMi,th) to the observed amplitude ratios 
(A?,obsA4i,obs)- The identified degree £ is the value 
which minimizes the % 2 : 



fe 

E 

J=2 L 



^4l,th ^4l,obs 

where k is the number of filters. 



(14) 



We note that the non-adiabatic predictions depend 
on some dominant parameters of the theoretical models 
(e.g. the metallicity for j3 Cephei stars and SPBs, the 
mixing length parameter a for S Scuti and 7 Doradus 
stars). Therefore, these parameters can be constrained 
by a feed-back process after a unique mode identification 
is achieved. We call this feed-back process non-adiabatic 
astero seismology , in which we iterate the procedure de- 
scribed above by adjusting the stellar parameters for the 
identified mode, until we find the best fit between theory 
and observations. We will illustrate this feed-back process 
below for the estimation of the metallicity of the (3 Cephei 
star EN Lac. 

4. Comparison with other methods 

The difference between our method and the one proposed 
by other authors is in the way of estimating the influ- 



ence of the effective temperature variation and the effec- 
tive gravity variation in Eqs. (j|) and ([b|). 

4.1. Mechanical boundary condition 

As a preliminary for the estimation of the relative effective 
temperature variation ST c g/T c f[ and the relative effective 
gravity variation Sg e /g ei the first step of the procedure fol- 



lowed by Wa tson ( |1988[ ), Heynderic kx et al. (g_994j), Cugier 
et al. ( |L994 ) and Balona & Evers ( 1999] ) was to compute 
SP/P at the photosphere. As initially proposed by Buta 
& Smith (1979), they used the following formula: 



SP_ 
~P 



— = [£(£ + l)K 



K ^ 



k 
Ro 



where K (sometimes denoted by ah) is given by 
GM 



K 



an 



(15) 



(16) 



Eq. (15) is deduced fr om th e "classical" mechanical 
boundary condition (Cox 1980 , Eq. (17.69')), by neglect- 



ing the Eulerian variation of the gravitational potential 
at the photosphere (Cowling approximation). The advan- 
tage of Eq. (15) is that it gives SP/P, without having to 
compute numerically the adiabatic or non-adiabatic eigen- 
functions throughout the entire star. Assuming that: 

^-(8,<j),t) = - CePp (cos 6) cos(at + m<jy) , (17) 



we have therefore: 



C = 4 



K~ 



£(£ + 1) K . 



(18) 



4.2. Influence of the effective temperature variation 



Wat son (|1988D , Garrido et al. ( |1990D and Heynderickx et 
al. (1994) related the Lagrangian temperature variation 



to the Lagrangian pressure variation at the photosphere 
by introduced a free ad-hoc parameter R describing the 
departure from adiabatic conditions. Moreover, they as- 
sumed that the Lagrangian temperature variation is equal 
to the local effective temperature variation at the photo- 
sphere: 



5T cS 



ST 



(19) 



The coefficient fx is then given by: 



r 2 - 1 

R^—\C\ 



= R 



r 2 - 1 



4 + R- 1 - £{£+l)K\ 



(20) 



In the adiabatic case, R = 1. Concerning the phase-lag 
ipT, Heynderickx et al. (1994) take the adiabatic value of 
180° for the applications to (3 Cephei stars and Garrido 
et al. (1990) let i/jt be a free parameter between 90° and 
135° for the applications to 5 Scuti stars. 
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Our approach does not make any of these assump- 
tions. We rigorously compute both the amplitude and the 
phase of the local effective temperature variation by non- 
adiabatic computations throughout the entire star and, in 
particular, throughout the entire non-grey atmosphere. 

Cugier et al. fll994j ) and Balona & Evers fll999| ) per- 
formed non-adiabatic computations in order to determine 
fx and tpT in a more rigorous way. In their method, they 
assume also that the Lagrangian temperature variation is 
equal to the local effective temperature variation at the 
photosphere (Eq. (19)). We do not make this assump- 
tion in our method, since we have shown in Dupret et 



al. ( |2002a ) that the Lagrangian variation of the tempera- 
ture at the photosphere is different from the variation of 
the local effective temperature, because of the significant 
Lagrangian variation of the optical depth. 

4.3. Influence of the effective gravity variation 



Stamf ord & Watson (|198l|) , Wat son (1988| ), Garrido et al. 
(|1990D and Heynderickx et al. ( |1994 ) related the coeffi- 
cient f g corresponding to the effective gravity variation 
(i.e. the gravity variation corrected for the pulsational ac- 
celeration) to the Lagrangian pressure variation by the 
following equation: 



fg = P* C = 



d In g 
dhiP 



C. 



(21) 



where p* is computed from equilibri um a tmosphere mod- 
els such as th e mod els of Kurucz ( 1993). Some authors 



(Cugier et al. [1994 
take p 



Balona & Evers 1999 ) proposed to 



1. 



Dupret et al. (2002a, Eq. (24)) proposed a more accu- 
rate way to determine Sg e /g e . In our method, f g is given 
by: 



Ro dtp'/dr 



M r 



K- 



(22) 



Under the Cowling approximation and neglecting the den- 
sity at the photosphere compared to the mean density of 
the star, Eq. (|22|) gives: 



fg = 2 



K~ 



(23) 



The correction leading to Eq. ( p3[ ) was also proposed 
by Cugier & Daszynska (2001). The difference between 
Eq. (|23| ) and Eq. (|2l| ) with p* = 1 is due to the fact that 
the Lagrangian variation of surface elements of the pho- 
tosphere (2 — £(£ + 1)K) affects the Lagrangian variation 
of the pressure described in Eq. (^8|), but does not affect 
the Lagrangian variation of the effective gravity. A simple 
comparison shows that: 

For p-modes : K is small so that the difference between 
Eq. ( plf ) and Eq. (^3|) is approximately 2. 

For high-order g-modes : K is large, so that the difference 
between Eq. (^) and Eq. ( p3| ) becomes very impor- 
tant ! 



Therefore, our improvement in the determination of the 
effective gravity variation has the largest impact on the 
photometric mode identification of g-mode pulsators such 
as SPBs and 7 Doradus stars. 

5. Applications 

In this section, we present the application of our mode 
identification method to two SPBs and one j3 Cephei star. 
The theoretical stellar models we used have been com- 
puted by the Code Liegeois d'Evolution Stellaire (CLES). 



5.1. Slowly pulsating B stars 

The two SPBs for which we performed non-adiabatic 
computations and a photometric mode identification are 
HD 74560 and HD 138764. Data obtained with Geneva 



photometry were taken from De Cat & Aerts (2002). We 
give in Table [l] the effective temperature and the gravity 
of these two stars, as derived from the most recent cali- 



bration of Kiinzli et al. (1997). We selected then theoret- 
ical models closest to these observations and their global 
characteristics are given in Tables || and ^. Numerous ad- 
ditional applications to other SPBs will be presented in 
De Cat et al. (in preparation). 

Table 1. Effective temperature and gravity of the stars 
HD 74560 and HD 138764 as deduced from Geneva cali- 
brations of Kiinzli et al. 



(|1997| ) 



HD 74560 Teft = 16210 ± 150 K logg = 4.15 ± 0.14 
HD 138764 T cff = 14050 ± 80 K \ogg = 4.20 ± 0.12 



Table 2. Global characteristics of the theoretical model 
of HD 74560 



M/M Q = 4.9 T cff = 16205 K log(L/L ) = 2.7521 
\ogg = 4.1677 R/Rq = 3.0208 age (My) = 24.5 
X = 0.7 Z = 0.02 no overshooting 



In Figs. and § we give the values of fx (local effective 
temperature variation for a normalized radial displace- 
ment at the photosphere) and in Figs. ^| and |4| those of ipx 
(phase difference between the local effective temperature 
variation and the radial displacement at the photosphere, 
in degrees) for different modes of the two SPBs, as a func- 
tion of the pulsation frequency in cd _1 , as computed by 
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Table 3. Global characteristics of the theoretical model 
of HD 138764 



M/Mq = 3.9 T efi = 14047 K log(L/L ) = 2.3760 
logg = 4.1964 R/R Q = 2.6073 age (My) = 38 
X = 0.7 Z = 0.02 no overshooting 



Table 4. Non-adiabatic results and mode identification 
for the star HD 74560. Degree £, radial order, amplitude 
of local effective temperature variation fx and phase-lag 
ipx for the modes with theoretical frequency closest to 
the observed dominant frequency / = 0.64472 cd _1 . The 
identified mode is given in bold. 



£ 


g n 


h 


lM°) 


1 


Sir 


10.08 


-23.0 


2 


g29 


22.38 


-8.2 


3 


g41 


36.01 


1.3 



Table 5. Non-adiabatic results and mode identification 
for the star HD 138764. Degree £, radial order, fx and 
tpx for the modes with theoretical frequency closest to 
the observed dominant frequency / = 0.7944 cd" 1 . The 
identified mode is given in bold. 



£ 


g n 


It 


It (°) 


1 


gl6 


5.59 


-23.2 


2 


g28 


17.23 


-9.3 


3 


g40 


31.57 


1.4 



our non-adiabatic code. In these figures, the vertical line 
corresponds to the observed frequency of the dominant 
mode. We see that, for a given frequency, the amplitude 
and phase-lag depend strongly on the degree I. The phys- 
ical explanation is the following. The term corresponding 
to the transversal compression in the equation of mass 
conservation is proportional to £(£ + 1). For high-order 
g-modes, this term dominates, which implies a strong de- 
pendence on £ of the eigenfunctions. The phase-lags (i/jt) 
are relatively close to zero, which is in agreement with the 
observations. 

In Tables | and § we give for the two stars, the degree 
£, the radial order, fx and ipx for the modes with the- 
oretical frequency closest to the observed frequency. The 
identified mode is given in bold. 
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Fig. 1. fx (local effective temperature variation for a nor- 
malized radial displacement at the photosphere) as func- 
tion of the pulsation frequency (in cd _1 ), for different 
modes of the SPB star HD 74560. The "+" correspond to 
modes of degree £ = 1, the "x" correspond to modes of 
degree 1 = 2 and the asterisks correspond to modes of de- 
gree I = 3. The vertical line corresponds to the observed 
frequency of the dominant mode. 
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Fig. 2. ipx (phase difference between the local effective 
temperature variation and the radial displacement at the 
photosphere in degrees) as function of the pulsation fre- 
quency in cd _1 , for different modes of the SPB star HD 
74560. The "+" correspond to modes of degree £ = 1, the 
" x " correspond to modes of degree £ = 2 and the aster- 
isks correspond to modes of degree £ = 3. The vertical 
line corresponds to the observed frequency of the domi- 
nant mode. 

In Figs. and ^, we give the amplitude ratios obtained 
from Geneva photometry for the dominant modes of the 
two stars. The bullets with error bars correspond to the 
observations. The full lines correspond to the theoretical 
predictions for different degrees £: solid line for £ = 1, 
dashed line for £ = 2 and dot-dashed line for I = 3. For 
both stars, a solution inside the error bars is found, and 
the identified degree is i = 1. These photometric mode 
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for different modes of the SPB star HD 138764. The "+" 
correspond to modes of degree I = 1, the " x " correspond 
to modes of degree 1 = 1 and the asterisks correspond to 
modes of degree t = 3. The vertical line corresponds to 
the observed frequency of the dominant mode. 



10 

5 

-5 
-10 
-15 
-20 
-25 



0.65 0.7 



0.75 



1.05 



0.8 0.85 0.9 0.95 

frequency (c/d) 

Fig. 4. ipT (in degrees) as function of the pulsation fre- 
quency in cd -1 , for different modes of the SPB star 
HD 138764. The "+" correspond to modes of degree I = 1, 
the " x " correspond to modes of degree £ = 2 and the as- 
terisks correspond to modes of degree £ = 3. The vertical 
line corresponds to the observed frequency of the domi- 
nant mode. 

identifications are in very good agreement with the spec- 
troscopic mode identifications performed by De Cat et al. 
(in preparation) using the moment method. 

5.2. The (3 Cephei star EN (16) Lac 

We present now the application of our method to the 
(3 Cephei star EN Lac. This star has been studied by 
many authors. We refer to Chapellier et al. ( 1995| ) and 
Lchmann et al. (2001) for a summary of the observational 
studies and to Dziembowski & Jerzykiewicz (1996) for the 
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Fig. 5. Amplitude ratios obtained with Geneva photom- 
etry for the dominant mode of the SPB star HD 74560. 
The bullets with error bars correspond to the observa- 
tions. The lines correspond to the theoretical predictions 
for different degrees £: solid line for £ = 1, dashed line for 
I — 2 and dot-dashed line for £ — 3. 
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Fig. 6. Amplitude ratios obtained with Geneva photom- 
etry for the dominant mode of the SPB star HD 138764, 
same caption as in Fig. ^| 

we illustrate the process we term non-adiabatic astero- 
seismology, by deriving constraints on the metallicity of 
the star. The three observed frequencies used in our study 



were taken from Lehmann et al. (2001) and the photomet- 
ric amplitudes obtained with Johnson filters were derived 



by Jerzykiewicz (1993) 



first seismic study of this (i Cephei star. In this section, 



Acrts et al. ( 2002 ) showed convincingly that the spec- 
troscopic mode identification is fully compatible with the 
photometric one for this star, and points towards a radial 
mode for the first frequency (/i = 5.9112 cd -1 ), an £ = 2 
mode for the second frequency (/2 = 5.8551 cd -1 ) and an 
I = 1 mode for the third frequency (/3 = 5.5033 cd -1 ). 
We therefore adopt this result here. We present the results 
of the amplitude ratios obtained for models with different 
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Table 6. Global characteristics of the theoretical models 
of EN Lac 
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metallicities. The choice of the models has been made with 
the following procedure. We computed models with 3 dif- 
ferent metallicities : Z = 0.015, Z = 0.02 and Z = 0.025. 
Because of the uncertainties in the calibration of the effec- 
tive temperature of EN Lac (see Jerzykiewicz & Sterken 



198C, Shobbrook 1985 and Dziembowski & Jerzykiewicz 
1996), we computed for each metallicity two evolutionary 
tracks with two different masses. For each evolutionary 
track, we subsequently selected the model giving the best 
agreement between the theoretical and observed frequen- 
cies, relying on the unambiguous mode identification. In 
all these models, X = 0.7 and there is no overshooting. 
Their global characteristics are given in Table [| 

In Fig. 0, we present the values of fr as a function of 
the pulsation frequency in cd -1 , for different modes and 
for the six models of EN Lac given in Table |[ The three 
vertical lines correspond to the three observed frequencies. 
We see that, the higher the metallicity, the lower the am- 
plitude of the local effective temperature variation for a 
normalized radial displacement (/r)- The physical origin 
of this phenomenon is explained in Fig. |^, where we show 
the amplitudes of the luminosity variation \8L/L\ as a 
function of the logarithm of temperature, from the center 
to the surface of the star, for the radial fundamental mode 
and for the models la, 2a and 3a of Table ^. The higher 
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Fig. 7. /t (local effective temperature variation for a nor- 
malized radial displacement at the photosphere) as func- 
tion of the pulsation frequency in c d _1 , for different modes 
(0 < i < 2) and for the six different models of the star 
EN Lac given in Table [| The "x" correspond to modes 
of the model la, the "+" correspond to the model lb, the 
full and empty squares correspond to the models 2a and 
2b respectively, the full and empty circles correspond to 
the models 3a and 3b respectively. The three vertical lines 
correspond to the three observed frequencies: v\ = 5.9112 
cd~\ v 2 = 5.8551 cd^ 1 and v 3 = 5.5033 cd" 1 . 



the metallicity, the more efficient the k mechanism, which 
implies a more important decrease of the luminosity varia- 
tion in the driving region. Therefore, the amplitude of the 
luminosity variation and of the local effective temperature 
variation at the photoshcre are smaller for a normalized 
displacement. The phase differences between the local ef- 
fective temperature variation and the radial displacement 
at the photosphere we obtained for the different models 
are very close to 180° (in agreement with the observations) 
and we do not give them here. By comparing in Fig. ^ the 
results obtained for the cold models (models la, 2a and 
3a) and the hot models (models lb, 2b and 3b), we see 
that changing the values of T e g within its observational 
error bars (keeping the metallicity constant) has only a 
very small effect on the non-adiabatic results. 

We present in Fig. ^ the theoretical amplitude ratios 
(Johnson photometry) obtained for three models of EN 
Lac with different metallicities (Table |6[ models la, 2a 
and 3a), and for the fundamental radial mode. We see 
that the model with Z = 0.015 gives the best agreement 
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Fig. 8. Amplitude of luminosity variation \8L/L\ as func- 
tion of the logarithm of temperature for the fundamental 
radial mode, for models of EN Lac with different metal- 
licities (models la, 2a and 3a of Table |J). The solid line 
corresponds to the model with Z = 0.015, the dashed line 
to the model with Z = 0.02 and the dotted line to the 
model with Z = 0.025. 
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between the theoretical and observed amplitude ratios. We 
have checked explicitly that all the modes in the observed 
range of frequencies remain unstable for this low metallic- 
ity; lower values are not compatible with mode excitation. 

The confrontation between the theoretical and ob- 
served amplitude ratios can thus be used as a constraint 
on the metallicity of stars driven by the metal opacity 
bump (/3 Cephei and Slowly Pulsating B stars), once we 
know the identification of the mode. We have seen in Fig. [?] 
that the non-adiabatic predictions and thus the theoreti- 



cal amplitude ratios are little affected by the uncertainties 
on T e ff for a given metallicity, so that the constraints we 
derived on the metallicity are reliable. This way of de- 
termining the metallicity may even turn out to be more 
precise than the classical method based on the analysis of 
the spectrum. We plan to validate our method to derive 
the metallicity by this feed-back process by confronting 
our predictions to those of j3 Cephei stars for which the 
metallicity is known with high accuracy. 

We note that, for rapidly rotating (3 Cephei stars, the 
interaction between pulsation and rotation can affect sig- 
nificantly the photometric amplitudes and phase-lags as 
shown recently by Daszynska et al. (2002). This interac- 
tion was not yet taken into account in our current non- 
adiabatic treatment. 



6. Conclusions 

We have presented an improvement of the often-used pho- 
tometric mode identification method. Our version of this 
method is based on precise non-adiabatic computations 
in which special attention is paid to the treatmen t of th e 
pulsation in the stellar atmosphere (Dupret et al. 2002a ). 
We have applied our new version of the method to identify 
the main mode of two SPBs (HD 74560 and HD 138764) 
and one (3 Cephei star (EN Lac, see also Aerts et al. 
|2002D . In both cases our photometric mode identifications 
were in very good agreement with the spectroscopic mode 
identifications, which are far less sensitive to temperature 
variations (De Ridder et al 



2002). We have shown also 



that the confrontation between the non-adiabatic theo- 
retical predictions and the observations can give interest- 
ing constraints on the models. We have used the term 
non-adiabatic astero seismology for this feed-back process. 
More precisely, for (3 Cephei stars and SPBs, the non- 
adiabatic predictions are very sensitive to the metallicity, 
so that this parameter can be constrained for these stars 
once definite mode identification is achieved. For S Scuti 
and 7 Doradus stars, the non-adiabatic predictions are 
very sensitive to the characteristics of the thin superficial 



convection la yer (Ba lona & Evers[1999[ Moya et al. [2002 , 
Dupret et al. |2002b ). Therefore, it is to be expected that 



a feed-back process similar to the one we presented will 
lead to a significant improvement of our understanding of 
this convection layer. 
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